Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome.

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

The Impact of Long-Term Hypoxia on the Antioxidant Defense System in the Siberian Frog Rana amurensis.

The Siberian frog Rana amurensis has a uniquely high tolerance to hypoxia among amphibians, as it is able to withstand several months underwater with almost no oxygen (0.2 mg/liter) vs. several days for other studied species. Since it was hypothesized that hypoxia actives the antioxidant defense system in hypoxia-tolerant animals, one would expect similar response in R. amurensis. Here, we studied the effect of hypoxia in the Siberian frog based on the transcriptomic data, activities of antioxidant enzyme, and content of low-molecular-weight antioxidants. Exposure to hypoxia upregulated expression of three relevant transcripts (catalase in the brain and two aldo-keto reductases in the liver). The activities of peroxidase in the blood and catalase in the liver were significantly increased, while the activity of glutathione S-transferase in the liver was reduced. The content of low-molecular-weight antioxidants (thiols and ascorbate) in the heart and liver was unaffected. In general, only a few components of the antioxidant defense system were affected by hypoxia, while most remained unchanged. Comparison to other hypoxia-tolerant species suggests species-specific adaptations to hypoxia-related ROS stress.

HiConfidence: a novel approach uncovering the biological signal in Hi-C data affected by technical biases.

The chromatin interaction assays, particularly Hi-C, enable detailed studies of genome architecture in multiple organisms and model systems, resulting in a deeper understanding of gene expression regulation mechanisms mediated by epigenetics. However, the analysis and interpretation of Hi-C data remain challenging due to technical biases, limiting direct comparisons of datasets obtained in different experiments and laboratories. As a result, removing biases from Hi-C-generated chromatin contact matrices is a critical data analysis step. Our novel approach, HiConfidence, eliminates biases from the Hi-C data by weighing chromatin contacts according to their consistency between replicates so that low-quality replicates do not substantially influence the result. The algorithm is effective for the analysis of global changes in chromatin structures such as compartments and topologically associating domains. We apply the HiConfidence approach to several Hi-C datasets with significant technical biases, that could not be analyzed effectively using existing methods, and obtain meaningful biological conclusions. In particular, HiConfidence aids in the study of how changes in histone acetylation pattern affect chromatin organization in Drosophila melanogaster S2 cells. The method is freely available at GitHub: https://github.com/victorykobets/HiConfidence.

SIRT6 is a key regulator of mitochondrial function in the brain.

The SIRT6 deacetylase has been implicated in DNA repair, telomere maintenance, glucose and lipid metabolism and, importantly, it has critical roles in the brain ranging from its development to neurodegeneration. Here, we combined transcriptomics and metabolomics approaches to characterize the functions of SIRT6 in mouse brains. Our analysis reveals that SIRT6 is a central regulator of mitochondrial activity in the brain. SIRT6 deficiency in the brain leads to mitochondrial deficiency with a global downregulation of mitochondria-related genes and pronounced changes in metabolite content. We suggest that SIRT6 affects mitochondrial functions through its interaction with the transcription factor YY1 that, together, regulate mitochondrial gene expression. Moreover, SIRT6 target genes include SIRT3 and SIRT4, which are significantly downregulated in SIRT6-deficient brains. Our results demonstrate that the lack of SIRT6 leads to decreased mitochondrial gene expression and metabolomic changes of TCA cycle byproducts, including increased ROS production, reduced mitochondrial number, and impaired membrane potential that can be partially rescued by restoring SIRT3 and SIRT4 levels. Importantly, the changes we observed in SIRT6-deficient brains are also occurring in aging human brains and particularly in patients with Alzheimer's, Parkinson's, Huntington's, and Amyotrophic lateral sclerosis disease. Overall, our results suggest that the reduced levels of SIRT6 in the aging brain and neurodegeneration initiate mitochondrial dysfunction by altering gene expression, ROS production, and mitochondrial decay.

Murine Falcor/LL35 lncRNA Contributes to Glucose and Lipid Metabolism In Vitro and In Vivo.

Glucose and lipid metabolism are crucial functional systems in eukaryotes. A large number of experimental studies both in animal models and humans have shown that long non-coding RNAs (lncRNAs) play an important role in glucose and lipid metabolism. Previously, human lncRNA DEANR1/linc00261 was described as a tumor suppressor that regulates a variety of biological processes such as cell proliferation, apoptosis, glucose metabolism and tumorigenesis. Here we report that murine lncRNA Falcor/LL35, a proposed functional analog of human DEANR1/linc00261, is predominantly expressed in murine normal hepatocytes and downregulated in HCC and after partial hepatectomy. The application of high-throughput approaches such as RNA-seq, LC-MS proteomics, lipidomics and metabolomics analysis allowed changes to be found in the transcriptome, proteome, lipidome and metabolome of hepatocytes after LL35 depletion. We revealed that LL35 is involved in the regulation of glycolysis and lipid biosynthesis in vitro and in vivo. Moreover, LL35 affects Notch and NF-κB signaling pathways in normal hepatocytes. All observed changes result in the decrease in the proliferation and migration of hepatocytes. We demonstrated similar phenotype changes between murine LL35 and human linc00261 depletion in vitro and in vivo that opens the opportunity to translate results for LL35 from a liver murine model to possible functions of human lncRNA linc00261.

Comparison of genome architecture at two stages of male germline cell differentiation in Drosophila.

Eukaryotic chromosomes are spatially segregated into topologically associating domains (TADs). Some TADs are attached to the nuclear lamina (NL) through lamina-associated domains (LADs). Here, we identified LADs and TADs at two stages of Drosophila spermatogenesis - in bamΔ86 mutant testes which is the commonly used model of spermatogonia (SpG) and in larval testes mainly filled with spermatocytes (SpCs). We found that initiation of SpC-specific transcription correlates with promoters' detachment from the NL and with local spatial insulation of adjacent regions. However, this insulation does not result in the partitioning of inactive TADs into sub-TADs. We also revealed an increased contact frequency between SpC-specific genes in SpCs implying their de novo gathering into transcription factories. In addition, we uncovered the specific X chromosome organization in the male germline. In SpG and SpCs, a single X chromosome is stronger associated with the NL than autosomes. Nevertheless, active chromatin regions in the X chromosome interact with each other more frequently than in autosomes. Moreover, despite the absence of dosage compensation complex in the male germline, randomly inserted SpG-specific reporter is expressed higher in the X chromosome than in autosomes, thus evidencing that non-canonical dosage compensation operates in SpG.

The Hitchhiker's Guide to Untargeted Lipidomics Analysis: Practical Guidelines.

Lipidomics is a newly emerged discipline involving the identification and quantification of thousands of lipids. As a part of the omics field, lipidomics has shown rapid growth both in the number of studies and in the size of lipidome datasets, thus, requiring specific and efficient data analysis approaches. This paper aims to provide guidelines for analyzing and interpreting lipidome data obtained using untargeted methods that rely on liquid chromatography coupled with mass spectrometry (LC-MS) to detect and measure the intensities of lipid compounds. We present a state-of-the-art untargeted LC-MS workflow for lipidomics, from study design to annotation of lipid features, focusing on practical, rather than theoretical, approaches for data analysis, and we outline possible applications of untargeted lipidomics for biological studies. We provide a detailed R notebook designed specifically for untargeted lipidome LC-MS data analysis, which is based on xcms software.

Long-Term Fluoxetine Administration Causes Substantial Lipidome Alteration of the Juvenile Macaque Brain.

Fluoxetine is an antidepressant commonly prescribed not only to adults but also to children for the treatment of depression, obsessive-compulsive disorder, and neurodevelopmental disorders. The adverse effects of the long-term treatment reported in some patients, especially in younger individuals, call for a detailed investigation of molecular alterations induced by fluoxetine treatment. Two-year fluoxetine administration to juvenile macaques revealed effects on impulsivity, sleep, social interaction, and peripheral metabolites. Here, we built upon this work by assessing residual effects of fluoxetine administration on the expression of genes and abundance of lipids and polar metabolites in the prelimbic cortex of 10 treated and 11 control macaques representing two monoamine oxidase A (MAOA) genotypes. Analysis of 8871 mRNA transcripts, 3608 lipids, and 1829 polar metabolites revealed substantial alterations of the brain lipid content, including significant abundance changes of 106 lipid features, accompanied by subtle changes in gene expression. Lipid alterations in the drug-treated animals were most evident for polyunsaturated fatty acids (PUFAs). A decrease in PUFAs levels was observed in all quantified lipid classes excluding sphingolipids, which do not usually contain PUFAs, suggesting systemic changes in fatty acid metabolism. Furthermore, the residual effect of the drug on lipid abundances was more pronounced in macaques carrying the MAOA-L genotype, mirroring reported behavioral effects of the treatment. We speculate that a decrease in PUFAs may be associated with adverse effects in depressive patients and could potentially account for the variation in individual response to fluoxetine in young people.

FACS-based isolation of fixed mouse neuronal nuclei for ATAC-seq and Hi-C.

The organization of chromatin structure plays a crucial role in gene expression, DNA replication, and repair. Chromatin alterations influence gene expression, and modifications could be associated with genomic instability in the cells during aging or diseases. Here, we provide a modified protocol to isolate fixed neuronal nuclei from a single mouse cortex to investigate the spatial organization of chromatin structure on a genome-wide scale by ATAC-seq (the assay for transposase-accessible chromatin with high-throughput sequencing) and chromatin conformation by Hi-C (high-throughput chromosome conformation capture).

IGF2 Mediates Resistance to Isoform-Selective-Inhibitors of the PI3K in HPV Positive Head and Neck Cancer.

Over 50% of human papilloma positive head-and-neck cancer (HNCHPV+) patients harbor genomic-alterations in PIK3CA, leading to hyperactivation of the phosphatidylinositol-4, 5-bisphosphate 3-kinase (PI3K) pathway. Nevertheless, despite PI3K pathway activation in HNCHPV+ tumors, the anti-tumor activities of PI3K pathway inhibitors are moderate, mostly due to the emergence of resistance. Thus, for potent and long-term tumor management, drugs blocking resistance mechanisms should be combined with PI3K inhibitors. Here, we delineate the molecular mechanisms of the acquisition of resistance to two isoform-selective inhibitors of PI3K (isiPI3K), alpelisib (BYL719) and taselisib (GDC0032), in HNCHPV+ cell lines. By comparing the transcriptional landscape of isiPI3K-sensitive tumor cells with that of their corresponding isiPI3K-acquired-resistant tumor cells, we found upregulation of insulin growth factor 2 (IGF2) in the resistant cells. Mechanistically, we show that upon isiPI3K treatment, isiPI3K-sensitive tumor cells upregulate the expression of IGF2 to induce cell proliferation via the activation of the IGF1 receptor (IGF1R). Stimulating tumor cells with recombinant IGF2 limited isiPI3K efficacy and released treated cells from S phase arrest. Knocking-down IGF2 with siRNA, or blocking IGF1R with AEW541, resulted in superior anti-tumor activity of isiPI3K in vitro and ex vivo. In vivo, the combination of isiPI3K and IGF1R inhibitor induced stable disease in mice bearing either tumors generated by the HNCHPV+ UM-SCC47 cell line or HPV+ patient-derived xenografts. These findings indicate that IGF2 and the IGF2/IGF1R pathway may constitute new targets for combination therapies to enhance the efficacy of PI3K inhibitors for the treatment of HNCHPV+.

TP73-AS1 is induced by YY1 during TMZ treatment and highly expressed in the aging brain.

Aging is a factor associated with poor prognosis in glioblastoma (GBM). It is therefore important to understand the molecular features of aging contributing to GBM morbidity. TP73-AS1 is a long noncoding RNA (lncRNA) over expressed in GBM tumors shown to promote resistance to the chemotherapeutic temozolomide (TMZ), and tumor aggressiveness. How the expression of TP73-AS1 is regulated is not known, nor is it known if its expression is associated with aging. By analyzing transcriptional data obtained from natural and pathological aging brain, we found that the expression of TP73-AS1 is high in pathological and naturally aging brains. YY1 physically associates with the promoter of TP73-AS1 and we found that along with TP73-AS1, YY1 is induced by TMZ. We found that the TP73-AS1 promoter is activated by TMZ, and by YY1 over expression. Using CRISPRi to deplete YY1, we found that YY1 promotes up regulation of TP73-AS1 and the activation of its promoter during TMZ treatment. In addition, we identified two putative YY1 binding sites within the TP73-AS1 promoter, and used mutagenesis to find that they are essential for TMZ mediated promoter activation. Together, our data positions YY1 as an important TP73-AS1 regulator, demonstrating that TP73-AS1 is expressed in the natural and pathological aging brain, including during neurodegeneration and cancer. Our findings advance our understanding of TP73-AS1 expression, bringing forth a new link between TMZ resistance and aging, both of which contribute to GBM morbidity.

Reduced purine biosynthesis in humans after their divergence from Neandertals.

We analyze the metabolomes of humans, chimpanzees, and macaques in muscle, kidney and three different regions of the brain. Although several compounds in amino acid metabolism occur at either higher or lower concentrations in humans than in the other primates, metabolites downstream of adenylosuccinate lyase, which catalyzes two reactions in purine synthesis, occur at lower concentrations in humans. This enzyme carries an amino acid substitution that is present in all humans today but absent in Neandertals. By introducing the modern human substitution into the genomes of mice, as well as the ancestral, Neandertal-like substitution into the genomes of human cells, we show that this amino acid substitution contributes to much or all of the reduction of de novo synthesis of purines in humans.

Aging and pathological aging signatures of the brain: through the focusing lens of SIRT6.

Brain-specific SIRT6-KO mice present increased DNA damage, learning impairments, and neurodegenerative phenotypes, placing SIRT6 as a key protein in preventing neurodegeneration. In the aging brain, SIRT6 levels/activity decline, which is accentuated in Alzheimer's patients. To understand SIRT6 roles in transcript pattern changes, we analyzed transcriptomes of young WT, old WT and young SIRT6-KO mice brains, and found changes in gene expression related to healthy and pathological aging. In addition, we traced these differences in human and mouse samples of Alzheimer's and Parkinson's diseases, healthy aging and calorie restriction (CR). Our results define four gene expression categories that change with age in a pathological or non-pathological manner, which are either reversed or not by CR. We found that each of these gene expression categories is associated with specific transcription factors, thus serving as potential candidates for their category-specific regulation. One of these candidates is YY1, which we found to act together with SIRT6 regulating specific processes. We thus argue that SIRT6 has a pivotal role in preventing age-related transcriptional changes in brains. Therefore, reduced SIRT6 activity may drive pathological age-related gene expression signatures in the brain.

Order and stochasticity in the folding of individual Drosophila genomes.

Mammalian and Drosophila genomes are partitioned into topologically associating domains (TADs). Although this partitioning has been reported to be functionally relevant, it is unclear whether TADs represent true physical units located at the same genomic positions in each cell nucleus or emerge as an average of numerous alternative chromatin folding patterns in a cell population. Here, we use a single-nucleus Hi-C technique to construct high-resolution Hi-C maps in individual Drosophila genomes. These maps demonstrate chromatin compartmentalization at the megabase scale and partitioning of the genome into non-hierarchical TADs at the scale of 100 kb, which closely resembles the TAD profile in the bulk in situ Hi-C data. Over 40% of TAD boundaries are conserved between individual nuclei and possess a high level of active epigenetic marks. Polymer simulations demonstrate that chromatin folding is best described by the random walk model within TADs and is most suitably approximated by a crumpled globule build of Gaussian blobs at longer distances. We observe prominent cell-to-cell variability in the long-range contacts between either active genome loci or between Polycomb-bound regions, suggesting an important contribution of stochastic processes to the formation of the Drosophila 3D genome.

Cumulative contact frequency of a chromatin region is an intrinsic property linked to its function.

Regulation of gene transcription is a complex process controlled by many factors, including the conformation of chromatin in the nucleus. Insights into chromatin conformation on both local and global scales can be provided by the Hi-C (high-throughput chromosomes conformation capture) method. One of the drawbacks of Hi-C analysis and interpretation is the presence of systematic biases, such as different accessibility to enzymes, amplification, and mappability of DNA regions, which all result in different visibility of the regions. Iterative correction (IC) is one of the most popular techniques developed for the elimination of these systematic biases. IC is based on the assumption that all chromatin regions have an equal number of observed contacts in Hi-C. In other words, the IC procedure is equalizing the experimental visibility approximated by the cumulative contact frequency (CCF) for all genomic regions. However, the differences in experimental visibility might be explained by biological factors such as chromatin openness, which is characteristic of distinct chromatin states. Here we show that CCF is positively correlated with active transcription. It is associated with compartment organization, since compartment A demonstrates higher CCF and gene expression levels than compartment B. Notably, this observation holds for a wide range of species, including human, mouse, and Drosophila. Moreover, we track the CCF state for syntenic blocks between human and mouse and conclude that active state assessed by CCF is an intrinsic property of the DNA region, which is independent of local genomic and epigenomic context. Our findings establish a missing link between Hi-C normalization procedures removing CCF from the data and poorly investigated and possibly relevant biological factors contributing to CCF.

Single-cell-resolution transcriptome map of human, chimpanzee, bonobo, and macaque brains.

Identification of gene expression traits unique to the human brain sheds light on the molecular mechanisms underlying human evolution. Here, we searched for uniquely human gene expression traits by analyzing 422 brain samples from humans, chimpanzees, bonobos, and macaques representing 33 anatomical regions, as well as 88,047 cell nuclei composing three of these regions. Among 33 regions, cerebral cortex areas, hypothalamus, and cerebellar gray and white matter evolved rapidly in humans. At the cellular level, astrocytes and oligodendrocyte progenitors displayed more differences in the human evolutionary lineage than the neurons. Comparison of the bulk tissue and single-nuclei sequencing revealed that conventional RNA sequencing did not detect up to two-thirds of cell-type-specific evolutionary differences.

TGF-Beta-Activated Cancer-Associated Fibroblasts Limit Cetuximab Efficacy in Preclinical Models of Head and Neck Cancer.

Most head and neck cancer (HNC) patients are resistant to cetuximab, an antibody against the epidermal growth factor receptor. Such therapy resistance is known to be mediated, in part, by stromal cells surrounding the tumor cells; however, the mechanisms underlying such a resistance phenotype remain unclear. To identify the mechanisms of cetuximab resistance in an unbiased manner, RNA-sequencing (RNA-seq) of HNC patient-derived xenografts (PDXs) was performed. Comparing the gene expression of HNC-PDXs before and after treatment with cetuximab indicated that the transforming growth factor-beta (TGF-beta) signaling pathway was upregulated in the stromal cells of PDXs that progressed on cetuximab treatment (CetuximabProg-PDX). However, in PDXs that were extremely sensitive to cetuximab (CetuximabSen-PDX), the TGF-beta pathway was downregulated in the stromal compartment. Histopathological analysis of PDXs showed that TGF-beta-activation was detected in cancer-associated fibroblasts (CAFs) of CetuximabProg-PDX. These TGF-beta-activated CAFs were sufficient to limit cetuximab efficacy in vitro and in vivo. Moreover, blocking the TGF-beta pathway using the SMAD3 inhibitor, SIS3, enhanced cetuximab efficacy and prevented the progression of CetuximabProg-PDX. Altogether, our findings indicate that TGF-beta-activated CAFs play a role in limiting cetuximab efficacy in HNC.

A machine learning framework for the prediction of chromatin folding in Drosophila using epigenetic features.

Technological advances have lead to the creation of large epigenetic datasets, including information about DNA binding proteins and DNA spatial structure. Hi-C experiments have revealed that chromosomes are subdivided into sets of self-interacting domains called Topologically Associating Domains (TADs). TADs are involved in the regulation of gene expression activity, but the mechanisms of their formation are not yet fully understood. Here, we focus on machine learning methods to characterize DNA folding patterns in Drosophila based on chromatin marks across three cell lines. We present linear regression models with four types of regularization, gradient boosting, and recurrent neural networks (RNN) as tools to study chromatin folding characteristics associated with TADs given epigenetic chromatin immunoprecipitation data. The bidirectional long short-term memory RNN architecture produced the best prediction scores and identified biologically relevant features. Distribution of protein Chriz (Chromator) and histone modification H3K4me3 were selected as the most informative features for the prediction of TADs characteristics. This approach may be adapted to any similar biological dataset of chromatin features across various cell lines and species. The code for the implemented pipeline, Hi-ChiP-ML, is publicly available: https://github.com/MichalRozenwald/Hi-ChIP-ML.

Differences in lipidome and metabolome organization of prefrontal cortex among human populations.

Human populations, despite their overwhelming similarity, contain some distinct phenotypic, genetic, epigenetic, and gene expression features. In this study, we explore population differences at yet another level of molecular phenotype: the abundance of non-polar and polar low molecular weight compounds, lipids and metabolites in the prefrontal cortical region of the brain. We assessed the abundance of 1,670 lipids and 258 metabolites in 146 Han Chinese, 97 Western European, and 60 African American individuals of varying ages, covering most of the lifespan. The statistical analysis and logistic regression models both demonstrated extensive lipid and metabolic divergence of the Han Chinese individuals from the other two populations. This divergence was age-dependent, peaking in young adults, and involved metabolites and lipids clustering in specific metabolic pathways.

Metabolome signature of autism in the human prefrontal cortex.

Autism spectrum disorder (ASD) is a common neurodevelopmental disorder with yet incompletely uncovered molecular determinants. Alterations in the abundance of low molecular weight compounds (metabolites) in ASD could add to our understanding of the disease. Indeed, such alterations take place in the urine, plasma and cerebellum of ASD individuals. In this work, we investigated mass-spectrometric signal intensities of 1,366 metabolites in the prefrontal cortex grey matter of 32 ASD and 40 control individuals. 15% of these metabolites showed significantly different intensities in ASD and clustered in 16 metabolic pathways. Of them, ten pathways were altered in urine and blood of ASD individuals (Fisher test, p < 0.05), opening an opportunity for the design of new diagnostic instruments. Furthermore, metabolic measurements conducted in 40 chimpanzees and 40 macaques showed an excess of metabolite intensity differences unique to humans, supporting the hypothesized disruption of evolutionary novel cortical mechanisms in ASD.